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SUMMARY 

This paper presents a deformation plasticity analysis of the tip region 
of a blunted crack in plane strain. The power hardening material is incompres- 
sible both elastically and plastically, in order to simulate behavior of a 
stress freezing material above critical temperature. The study represents a 
full field, finite difference solution to the Mode I problem. Stress and dis- 
placement fields surrounding the crack tip are presented. The results of this 
study indicate that the maximum stress seen at the crack tip is Indeed limited 
and is determined by the tensile properties; however, the scale over which the 
stresses act is dependent on the loading. Comparisons are good between the 
forward crack tip displacement and micro-fractographic measurements of 
"stretch" zones observed in plane strain fracture toughness tests. 


INTRODUCTION 

In recent years Cherepanov (ref. 1), Rice (ref. 2,3), Hutchinson (ref. 4, 
5), and Rice and Rosengren (ref. 6) have shown the asymptotic behavior of 
stress and strain fields surrounding sharp crack tips in plane strain. Using 
these studies as a guide, full field solutions with finite elements have been 
obtained by Levy, Marcal, Ostergren and Rice (ref. 7) and Hilton and Hutchinson 
(ref. 8). These two studies give accurate near and far field behavior due to 
the inclusion of singular elements reflecting plasticity at the crack tip. 

Other numerical solutions by Marcal and King (ref, 9), Mendelson (ref. 10), 
Swedlow and coworkers (ref. 11,12) and Tuba (ref. 13) show qualitative features 
of the near field, but may not yield accurate stress field definition due to 
the large gradients there. 

In order to have an accurate description of the near field surrounding 
crack tips, and hence a good understanding of the mechanisms of failure, 

Rice and Johnson (ref. 14) have pointed out that crack tip blunting must also 
be included. Their analysis accounted in an approximate manner for the 
blunting at the crack tip and for strain hardening in the plastic zone. As a 
result they showed that the stresses near the crack tip are Indeed finite and 
that the maximum Oyy stress occurred at some small distance from the deformed 
track tip. A finite deformation analysis by McGowan and Smith (ref. 15) of 
jlunted cracks in a linear (stress-strain) incompressible material shows the 
same general behavior. The maximum Oyy stress occurs in front of the blunted 
track tip and the magnitude is independent of the remote loading. 
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The purpose of the present study is to gain a full field solution around 
a blunted crack tip in a strain hardening incompressible material under Mode I 
loading. This work will provide an accurate description of the stress and de- 
formation fields immediately surrounding the blunted tip, and thereby gain 
insight to fracture behavior. Deformation theory of plasticity with a Mlses 
yield condition is used. The resulting set of equations is solved for the 
blunted crack tip in the deformed state under load by finite differences. 

The linear theory of Inglis (ref. 16) gives the necessary asymptotic boundary 
conditions . 

An initial goal of the present study was to gain a more complete under- 
standing of the near field behavior of stress freezing photoelastic materials 
above critical temperature; however, this study should also give considerable 
insight to the general behavior of engineering materials under Mode I loading. 

SYMBOLS 


c One-half crack length 

E Young * s Modulus 

K Stress intensity factor 

n Strain hardening exponent 

r,9 Cylindrical coordinates measured from 
crack tip 

T Secant modulus = a /e„ 

e T 

u^ Displacement vector 

U Strain energy density 

X Distance in front of deformed crack 

tip 

Y Distance perpendicular to deformed 
crack tip 


Initial yield strain 
Effective plastic strain 

e_ Effective total strain = 

■ e + a /E 

P e 

p Deformed crack root radiui 


V Poisson's ratio 


a . . Stress tensor 

ij 

a Hoop stress 

yy 

a Effective stress 

' e '• 

Oq Tensile yield stress 


f Airy stress function 


e.. Strain tensor 
Ij 


a Constant in eq. (2) 

FORMULATION OF THE PROBLEM 


Using small strain deformation theory of plasticity for an incompressible 
(v= 1/2) material the governing equation for the field can be shown to be; 
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(The details of this analysis are given in ref. 17) 
For this study a Ramberg-Osgood material will be used 


f" ba[(a^/a^)^/“ - 1] 


or f = 1 + - a^/a^] 


( 2 ) 


where b = 0 if a < o 
e 0 

b = 1 if a > a 
e o 

Thus the governing equation (1) will be solved subject to the constitutive 
laws (eq. (2)). 


The geometry of the blunted cracks in the deformed state under Mode I 
loading will resemble small elliptical perforations as shown in figure 1. 
The size of the deformed crack tip root radius will be determined through 
Integration of the strain displacement relationships 


+ u. 


2e, 


i.j j,i ""Ij 

The affected strain hardening region will be divided into a small grid 
utilizing elliptical coordinates and the governing set of equations will be 
solved through the method of finite differences. At some distance from the 
deformed crack tip the linear solution of Inglis (ref. 16) will apply. The 
stress at the outer boundary of the inner strain hardening region will be then 
matched to the Inglis solution. The outer boundary will be enlarged until 
there is no change in the inner stress field. (A detailed description of the 
solution procedure is given in ref. 18.) 


PRESENTATION OF RESULTS 


The stress and displacement fields in the field surrounding a deformed 
crack tip in a strain hardening material which is Incompressible in both the 
elastic and plastic regions are examined. Strain hardening exponents of 0.2 
through 0.01 are presented. The range of initial yield strain values is from 
0.01 through 0.0001. The value of a in the effective stress-effective strain 
relationship, equation (2), is taken to be 1.0 in this study. (The authors 
have found that small changes in a and V do not influence the solution signifi- 
cantly.) The "linear" results reported here are those of Inglis (ref. 16) for 
a deformed crack tip in a linear material. The "singular" results are those 
corresponding to a crack which has no root radius in a linear material. 

The plastic zone shape for the smallest ellipse investigated (p = 0.0018 
(K/0^)2 ) is shown in figure 2. Note that with decreasing hardening (n 0) 
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the plastic zone grows in maximum extent and leans progressively in the direc~ 
tion of crack propagation. For comparison, the singular plastic zone from 
McClintock and Irwin (ref. 19) and that of Levy et al (ref. 7) for a non- 
hardening (n = 0) material are also shown. As shown by figure 2, the plastic 
zone shape predicted by McClintock and Irwin (ref. 19) is approached by the 
present study as n ^ «>. The difference between the plastic zone shape pre- 
dicted by Levy et al (ref. 7) and the present study for n = 0.01 is primarily 
due to the inclusion of blunting effects and use of v of 0.5 in the latter; 

the difference should be negligible as e -^0. 

o 

The effective stress a is shown versus the distance ahead of the de- 
formed crack tip in figure 3? This figure shows that the effective stress 
varies as (rii/n+l)“l in the plastic zone ahead of the tip. (The behavior for 
Other values of 6 is similar) . It can be shown that the strain energy has the 

(1+n) /n 

- 1 for a power hardening material 

o L - 

Therefore, the strain energy varies approximately as 1/r in the plastic zone. 
This was a key assumption in the analysis of Rice and Rosengren (ref . 6) and 
Hutchinson (ref. 4) . 

The Oyy stress in front of the crack tip is shown for various values of 
yield strain for n = 0.01 in figure 4. As shown in this figure, this stress 
is substantially reduced near the crack tip because of blunting and strain 
hardening, with the maximum value developed at some distance forward of the 
crack tip. (The stress distributions for other values of n is quite 
analogous.) The analysis of Rice and Johnson (ref. 14) gives the same quali- 
tative behavior; the correlation is believed to be quite reasonable in view 
of the several approximations involved. For a non-hardening material Rice 
(ref. 2) has shown that the maximum Oyy stress is 2.97 o^. This stress, as 
predicted by Levy et al (ref. 7), approaches this limit at the crack tip as 
shown in figure 4. The Oyy stress distribution of the present study in this 
figure reflects the presence of blunting and should coincide with the work of 
Levy et al (ref. 7) as -> 0. 

Figure 5 shows the variation of maximum Oyy stress with initial yield 
strain for varying hardening. As shown in the figure, blunting alone (the 
"linear" curve) forces the peak Oyy stress to be finite and the inclusion of 
finite deformations (ref. 15) reduces the magnitude somewhat. However, the 
effects of blunting and plasticity taken together are significant: the peak 

CTyy stress is reduced by a factor of 10 from that with blunting alone. From 
figure 5, one observes the peak Oyy stress to be So to lo depending upon n 
and cTq/E. The peak o y stress increases with n and decreases with Oq/E. (The 
large value of peak ©yy stress compared to the uniaxial yield stress, Oq, is 
believed due to the presence of triaxiality in the crack tip region.) 

The crack tip displacement in the direction of propagation (which is also 
the deformed crack root radius, p) is shown in figure 6 for varying initial 
yield strain and hardening exponent. The present study predicts that the for- 
ward crack tip displacement Increases with Oq/E and decreases with n. For 
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eomparison one-half the crack tip opening displacement calculated by Levy (ref. 
7) is shown. The forward crack tip displacement as predicted by the present 
study and the work of Levy et al (ref. 7) show parallel behavior, although they 
are separated by some distance. This disagreement is believed due to the shape 
of the crack tip being elliptical in the present study instead of cylindrical. 

Included also on figure 6 is the width of the "transition" or "stretch" 
zone which exists on the fracture surface between the cracked and the overload 
regions in fatigue. As Broek (ref. 20) has discussed, the depth of this tran- 
sition zone is the crack tip opening displacement, and, therefore the width is 
the tip forward displacement. 

Examination of the figure shows the correlation between the forward tip 
displacement and failure. The measurements of the stretch zone fall close to 
n = 0.2. For the steels and alumlntans shown values of n around 0.05 have been 
reported in references 20, 21 and 22. However, it is known that for this class 
of materials the value of n varies with plastic strain (ref. 23). For large 
plastic strain (e_ > 10%), the strain hardening exponent is close to 0.2 as 
shown by Jones and Brown (ref. 24) for 4340 steel. The strains in the tip 
region are clearly greater than 10% so that the agreement between the measure- 
ments and the analysis appears quite reasonable. The scatter band shown on 
the figure is an indication of the span of actual measurements Cauthors typi- 
cally report a 40% variation) . 


DISCUSSION 

Previously McGowan and Smith (ref. 15) performed a finite deformation 
analysis of the region surrounding deformed crack tips for a linear (stress- 
strain) material. The results of the finite deformation work showed that the 
maximiim Oyy stress occurred in front of the deformed crack tip. It was deter- 
mined that the stress distribution around the crack tip was "similar", in the 
sense that one stress distribution could be used to describe the response of 
the material under load. The size of the affected zone would depend upon 
the load and crack length through K. The self-similarity of the stress field 
was a direct result of the blunting process, and would be expected to remain 
as long as the affected zone stayed small with respect to the crack length, 
thickness, or any other in-plane dimension. 

* 

The behavior is quite similar for a power hardening material. The stress 
field is self-similar with the size of the affected zone varying with K. The 
maximum <jyy stress will only then be a function of the material properties 
E, n, and Gq. The stretching of the similar stress distribution will depend 
upon K as well as the other material properties. One may conjecture that 
failure would depend upon the growth in size of a critical dimension, such as 
plastic zone size, which increases with K, 

Wells (ref. 25) and others have used the crack opening displacement as a 
fracture criterion. Broek (ref. 20) has used this concept to correlate the 
depth of transition zones in aluminum with fracture toughness. The present 
study shows good correlation of fracture toughness and transition zone width. 
Krafft (ref. 26), Hahn and Rosenfleld (ref. 27) and Rice and Johnson Ctef. 14) 


589 



have all shown good correlation of plane strain fracture toughness with some 
minute particle size or process zone size for specific cases. 

SUMMARY AND RECOMMENDATIONS 

Following the pioneering studies of Hutchinson (ref . 4 ) , Rice and Rosen- 
gren (ref. 6), Levy et al (ref. 7) and Hilton and Hutchinson (ref. 8), the 
authors have obtained a full field deformation plasticity finite difference 
solution to the Mode I plane strain problem including the effects of blunting. 
The material was Incompressible in both the elastic and plastic regions, and 
followed a power hardening rule. Stress and displacement fields surrounding 
the deformed crack tip are presented, and are found to compare favorably both 
with the analysis of other investigators as well as experimental results. Be- 
cause of the Improved accuracy expected from a full field solution, it would 
be appropriate to Incorporate such a solution into theories concerning void 
coalescence and final instability. Efforts are currently being devoted to 
such an approach. 

REFERENCES 

1. Cherepanov, G. P,, "Crack Propagation in Continuous Media," Appl. Math . 
Mech . (PMM), Vol. 31, p. 476, 1967. 

2. Rice, J. R. , "A Path Independent Integral and the Approximate Analysis of 
Strain Concentrations by Notches and Cracks", J. Appl. Mech ., Vol. 35, p. 
379, 1968. 

3. Rice, J., R. , "Mathematical Analysis in the Mechanics of Fracture", Ch. 3 
of Fracture, An Advanced Treatise, Vol. II (H. Liebowitz, ed.). Academic 
Press, New York, 1968. 

4. Hutchinson, J. W. , "Singular Behavior at the End of a Tensile Crack in a 
Hardening Material", J. Mech. Phys. Solids , Vol. 16, p. 13, 1968. 

5. Hutchinson, J. W. , "Plastic Stress and Strain Fields at a Crack Tip", £. 
Mech. Phys. Solids , Vol, 16, p. 337, 1968. 

6. Rice, J. R. and Rosengren, G. F., "Plane Strain Deformation Near a Crack 
Tip in a Power Law Hardening Material", J. Mech. Phys. Solids , Vol. 16, 
p . 1 , 1968 . 

7. Levy, N., Marcal, P. V., Ostergren, W. J., and Rice, J. R. , "Small Scale 
Yielding Near a Crack in Plane Strain: A Finite Element Analysis", 

Int . J . Frac . Mech . , Vol. 7, No. 2, p. 143, 1971. 

8. Hilton, P. D. and Hutchinson, J, W. , "Plasticity Intensity Factors for 
Cracked Plates", Eng . Frac . Mech . , Vol, 3, p. 435, 1971. 

9. Marcal, P.V. and King, I.P., "Elastic-Plastic Analysis of Two-Dimensional 
Stress Systems By the Finite Element Method", Int . J . Mech Sci . , Vol. 9, 
p. 143, 1967. 

10. Mendelson, A., Plasticity; Theory and Application , The Macmillan Company, 
New York, 1968. 


590 



11. 


Swedlow, J. L. , Yang, A. H., and Williams, M. L., "Elasto-Plastic Stresses 
and Strains in a Cracked Plate”, in Proc. of the First Int. Conf. on Frac- 
ture , 1965. Vol. 1, p. 259, 1966. 

12. Swedlow, J. L., "Elasto-Plastic Cracked Plates in Plane Strain”, Int. J . 
Frac. Mech ., Vol. 5, p. 33, 1969. 

13. Tuba, I. S., ”A Method of Elastic Plastic Plane Stress and Strain Analy- 
sis”, J. of Strain Analysis , Vbl. 1, p. 115, 1966. 

14. Rice, J. R. and Johnson, M. A., "The Role of Large Crack Tip Geometry 
Changes in Plane Strain Fracture”, Inelastic Behavior of Solids (M. F. 
Kanninen, et al, eds.), McGraw-Hill, New York, p. 641, 1970. 

15. McGowan, J. J. and Smith, C. W. , "A Finite Deformation Analysis of the 
Near Field Surrounding the Tip of Crack-Like Elliptical Perforations”, 

Int. J. Frac .. Vol. 11, No. 6, p. 977, 1975. 

16. Inglls, C. E., "Stresses in a Plate Due to the Presence of Cracks and 
Corners”, Trans. Instn. Naval Archit ., 55, p. 219, 1913. 

17. McGowan, J. J. and Smith, C. W. , "A Deformation Plasticity Analysis of the 
Blunted Crack Tip in Plane Strain,” VPI-E-76-4, March 1976. 

18. McGowan, J. J. and Smith, C. W. , "A Finite Deformation Analysis of the Near 
Field Surrounding the Tip of Crack-Like Ellipses”, VPI-E-74-10, May 1974. 

19. McClintock, F. A. and Irwin, G. R. , "Plasticity Aspects of Fracture Mech- 
anics”, Fracture Toughness Testing and Its Applications , STP-381, ASTM 

p. 84, 1965. 

20. Broek, D. , "Correlation Between Stretched Zone Size and Fracture Tough- 
ness”, Eng. Frac. Mech ., Vol. 6, No. 1, p. 173, 1974. 

21. Bates, R. C., Clark, Jr., W. G., and Moon, D. M., "Correlation of Fracto- 
graphic Features with Fracture Mechanics", Electron Microfractography, 
STP-453, ASTM, p. 192, 1969. 

22. Pandey, R. K. and Banerjee, S., "Studies on Fracture Toughness and Fracto- 
graphic Features in Fe-Jto Base Alloys", Eng. Frac. Mech ., Vol. 5, No. 4, 
p. 965, 1973. 

23. Lauta, F. J. and Steigerwald, E. A., "Influence of Work Hardening Coeffi- 
cient on Crack Propagation in High Strength Steels”, AFML TR 65-31, Air 
Force Materials Laboratory, May 1965. 

24. Jones, M. H. and Brown, Jr., W. B., "The Influence of Crack Length and 
Thickness in Plane Strain Fracture Toughness Tests", Review of Developments 
in Plane Strain Fracture Toughness Testing, STP-463, ASTM, p. 63, 1970. 

25. Wells, A. A., "Crack Opening Displacements from Elastic-Plastic Analyses of 
Externally Notched Tension Bars*', Eng. Frac. Mech ., VoL 1, No. 3, p. 3 99,1969. 

26. Krafft, J., "Correlation of Plane Strain Crack Toughness with Strain 
Hardening Characteristics of a Low, a Medium, and a High Strength Steel", 
Applied Material Research , Vol. 3, p. 1964, 1964. 

27. Hahn, G. and Rosenfield, A., "Source of Fracture Toughness: The Relation 

Between Kic and the Ordinary Tensile Properties of Metals", Applications 
Related Phenomena in Titaniinn Alloys, STP-432, ASTM, p. 5, 1968. 


591 






Figure 3.*- Effective stress distribution forward of the 
blunted crack tip. 



Figure hr Distribution of <Tyy stress forward of the blunted 
crack tip for n = 6.01. 







Forward crock lip displacement, p/{K/a^) 


Initial yield strain, S q 


Figure S.” Variation of CTyy atress maximum with and n. 



Figure 6.- Forward crack tip displacement and stretch zone 
variation with n and e . 

















